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Abstract 



We propose a microscopic effective interaction to treat pairing correlations in the 
^Sq cliannel. It is introduced by recasting tfie gap equation written in terms of the 
^ ' bare force into a fully equivalent pairing problem. Within this approach, the pro- 

2 ' posed interaction reproduces the pairing properties provided by the realistic Al^lS 

force very accurately. Written in the canonical basis of the actual Bogolyubov trans- 
formation, the force takes the form of an off-shell in-medium two-body matrix in the 
superfluid phase multiplied by a BCS occupation number 2 pm- This interaction is 
finite ranged, non local, total-momentum dependent and density dependent. The fac- 
tor 2 pm emerging from the recast of the gap equation provides a natural cut-off and 
makes zero-range approximations of the effective vertex meaningful. Performing such 
an approximation, the roles of the range and of the density dependence of the inter- 
action can be disentangled. The isoscalar and isovector density-dependences derived 
ab-initio provide the pairing force with a strong predictive power when extrapolated 
toward the drip-lines. Although finite ranged and non local, the proposed interaction 
makes HFB calculations of finite nuclei in coordinate space tractable. Through the 
two-basis method, its computational cost is of the same order as for a zero-range force. 
PACS: 21.60.-n; 21.30.-x; 
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1 Introduction 



The structure of the nucleus and the properties of extended nuclear systems strongly 
depend on their possible superfiuid nature. In finite nuclei, nucleonic pairing has a strong 
influence on all low-energy properties of the system. This encompasses masses, separation 
energies, deformation, individual excitation spectra and collective excitation modes such 
as rotation or vibration. The role of pairing correlations is particularly emphasized when 
going toward the drip-lines. This is due to the proximity of the Fermi surface to the 
single-particle continuum. Indeed, the scattering of virtual pairs into the continuum gives 
rise to a variety of new phenomena as regards the properties of ground and excited states 
of nuclei [1]. 

In neutron stars, a good description of pairing is also required. The neutron-neutron 
^5*0 pairing drives the physics of neutron star crusts, while at higher densities, that is from 
the inner crust to two or three times the saturation density, neutron and proton pairing 
occur dominantly in the ^i-*2 — '^-^2 and ^5*0 channels, respectively [2]. Such superfiuid 
phases infiuence the dynamical and thermal evolution of the star. Indeed, post-glitching 
timing observations [3] and the cooling history [4] strongly depend on the presence or 
absence of pairing in the system. 

To treat pairing, one needs to specify the many-body technique used and the appro- 
priate interaction to insert into the calculation at that chosen level of approximation. 
The latter depends on the situation and on the system. In the present case, our aim 
is eventually to perform non-relativistic self-consistent mean-field and beyond-mean-field 
calculations in finite nuclei. Mean-field calculations are of the Hartree-Fock-Bogolyubov 
(HFB) type, while the considered beyond-mean-field calculations deal typically with sym- 
metry restorations (Projected Mean Field Method) and with large amplitude motion in 
nuclei (Generator Coordinate Method) [5, 6]. Thus, one has to identify the appropriate 
vertices to be used in each of these cases. The same question arises for instance in the 
context of the shell model [7]. 

We concentrate here on the mean- field treatment. While the variational derivation of 
the HFB equations cannot help in defining the appropriate vertex, the Green function or 
Goldstone formalisms are able to do so. Such many-body theories show unambiguously 
that the irreducible vertex to be used in the pairing channel at lowest order is the bare 
nucleon-nucleon (NN) force [8, 9, 10, 11]. At the next order, the irreducible pairing vertex 
involves the so-called polarization diagrams. This situation is in contrast to the particle- 
hole channel where one cannot avoid regularizing the repulsive core of the bare interaction 
from the outset through the definition of an in-medium two-body matrix. This stresses 
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the fact that the effective forces in the particle-hole and particle-particle channels are 
different. In particular, the direct use in the gap equation of an in-medium vertex such as 
the Brueckner G-matrix leads to double counting [11, 12]. 

Thus, the mean-field energy defined in this context is a functional of one-body normal 
and abnormal density matrices and do not refer to the mean value of a given Hamiltonian 
in a product state. However, the strategy used consists of motivating the low-energy 
functional from a many-body expansion and thus, consists of keeping an explicit link to the 
bare NN force. With this approach, the theory is not necessarily local and the mean-field 
functional does not include more correlations than provided by the irreducible vertices in 
the particle-hole and particle-particle channels. Also, one has to come back to the many- 
body expansion to enrich the functional if going beyond the mean-field approach [13]. 
Even if the final goal is the same, this differs from the strategy used in the so-called (local) 
density functional theory [14, 15]. Note also that calculations based on the extension of 
Gorkov's formalism to the relativistic case are also performed [16, 17, 18]. Recently, the 
effects of polarization leading to the screening of nucleon and meson propagators have 
been studied in infinite matter [19]. The reduction of the gap was found to be much larger 
than in the non-relativistic case, correcting to some extent for the excessive gap found at 
first order. 

As regards HFB calculations in finite nuclei, only phenomenological forces have been 
used in the spin-singlet /isospin-triplet pairing channel so far. One example is the finite- 
range, density-independent Gogny force [20], whose restriction to the spin-singlet /isospin 
triplet channel can be written as: 

1 - P ^ 

yGosnj/^^^^^^) ^ ^ g-|r-i-f2|V«f , (1) 

1=1 

and is to be averaged over the angle between the incoming and outgoing relative momenta 
if dealing with the S wave only. In Eq. 1, is the spin-exchange operator while the 
Coulomb part of the Gogny force has not been considered. The other commonly used 
pairing interaction is the (density-dependent) delta interaction (DDDI) [21, 22, 23, 24, 25]: 

where p denotes the matter density (= local scalar-isoscalar part of the one-body density- 
matrix). 

The latter, usually used when solving the problem in coordinate space, must be com- 
plemented with a cut-off in the gap equation to avoid divergences. Studies of the rotational 
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bands of supcr-dcformcd nuclei [26, 22] and actinidcs [23, 27], of halo nuclei [21], and of 
the evolution of charge radii across magic numbers [28, 29] have helped to establish the 
success and the surface-peaked character of the DDDI in the pairing channel. Recently, 
more systematic studies of asymptotic matter and pair densities of exotic nuclei [30] , of the 
evolution of the pairing gap toward the neutron drip-line [31], and the average behavior 
of the odd-even mass differences over the mass table [32] have allowed a refinement of the 
DDDI. The optimal compatibility between experimental data and mean-field calculations 
was obtained for a force between surface and volume [23, 32], with pc ~ 2psat^ where psat = 
0.16 fm~^ is the saturation density of symmetric nuclear matter. The great sensitivity of 
matter and pair densities to pairing in the low-density regime seemed to favor 1/2 < 7 < 1. 

Although successful in describing low-energy nuclear structure over the (known) mass 
table [6, 22, 23, 33, 34, 35, 36, 37], the two previous phenomenological pairing interactions 
lack a link to the bare nucleon-nucleon interaction. They were directly fitted to finite nu- 
clei data, and may thus renormalize beyond-mean-field effects. In addition, their fits were 
performed where experimental data are available. Extrapolating the use of these interac- 
tions toward the drip- lines is questionable. To correct for this defect, few DDDI were fitted 
to reproduce the gap provided by realistic NN forces in infinite matter [38, 39]. However, 
the necessary density dependence and cut-off were still treated phenomenologically. 

In fact, the present knowledge about the pairing force and the nature of pairing cor- 
relations in nuclei is quite poor. Properties such as the range of the effective pairing 
interaction, its link to the bare force, its possible surface character in finite nuclei and its 
density dependence, in particular isovector, still have to be clarified. As noticed as early 
as thirty years ago [40] and pointed out several times since [1, 13, 29, 31, 32, 41, 42, 43], 
obtaining proper density dependences of particle-hole and particle-particle effective forces 
at a given level of approximation is difficult but of great importance to meet modern high- 
precision experiments. The problem dealing with the cut-off to be used in the pairing 
channel in connection with zero-range vertices has been solved recently [44, 45]. The idea 
was to identify the divergences stemming from the use of a local gap and regularize them 
through a well-defined renormalization scheme. That scheme has to be understood in the 
context of the local density functional theory [29] (which is perfectly fine) rather than as 
a mean-field approximation arising at lowest order of some many-body expansion. 

The present work concentrates on the mean-field treatment of pairing at relatively 
low energy and low density in the isotropic, spin-singlet and isospin-triplet channel. The 
aim of this study is manifold. In section 2, we define an appropriate simple version of 
the bare force in the channel and explain in details its fitting procedure. Then, 
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an in- medium pairing interaction (Eq. 21) equivalent to the bare force is introduced in 
section 3.1. The corresponding diagrammatic resummation authorizes the study of the 
finite-ranged effective force and of its zero-range approximation on the same footing. This 
is discussed in section 3.2 where the roles of the range and of the density dependence of the 
pairing interaction are disentangled. In the particular case of the zero-range approximation 
for the vertex, the scheme proposed presents strong similarities with the regularization 
procedure introduced in Refs. [44, 45]. We will actually discuss this particular point in a 
forthcoming publication. Ultimately, the interaction is to be used in calculations performed 
in coordinate space by solving the HFB equations on a three dimensional mesh [46, 47, 48]. 
Although finite ranged and non local, the proposed interaction is shown in section 3.3 to 
make these calculations tractable. The formulas defining completely the new effective 
pairing force can be found in the same section. Some important points are discussed in 
section 4 while our conclusions are given in section 5. 

2 Simplified bare force 
2.1 Fitting procedure 

Screening effects beyond the mean-field approximation due to density and spin fluctuations 
are known to strongly decrease the pairing gap in neutron matter, both for singlet and 
triplet pairing [2, 49, 50, 51, 52]. Whether it is justified to extend this statement to 
finite nuclei is still an open question. Indeed, the dressing of the vertex could change 
significantly when going from neutron matter to symmetric matter [53]. Also, including 
the induced interaction and off-shell self-energy effects due to the exchange of surface 
vibrations between time-reversed states seems to increase the pairing gap in finite nuclei 
compared to that generated by the bare force [26, 54, 55]. In addition, the influence of the 
restoration of particle-number and pairing vibrations still have to be characterized from 
systematic self-consistent calculations in even and odd nuclei. The situation as regards 
beyond-mean-field effects is unclear at this stage and one can simply state that a significant 
cancellation between the effects of screening and of surface vibrations on singlet pairing 
should take place in finite nuclei. As a result, two strategies seem reasonable when dealing 
with the definition of a pairing interaction to be used in mean-field calculations. 

Sticking to the pure mean-field picture, one can define the interaction by reproducing 
properties of the bare NN interaction. If necessary, the possibility remains to include 
beyond-mean-field effects explicitly in a consistent way when using an interaction mim- 
icking the bare force. 



5 



A second strategy eonsists of fitting tlic interaction directly to finite nuclei data through 
mean- field calculations [22]. This strategy has been the most popular so far when dealing 
with phenomenological forces to be used in self-consistent Hartree-Fock-Bardeen-Cooper- 
Schrieffer (HFBCS) and HFB calculations [20, 22]. Such a procedure aims at renormalizing 
the beyond-mean-field effects which possibly do not cancel out. Of course, the use of such 
a force in calculations going explicitly beyond the mean-field is suspicious. 

As a first attempt, and because we want to separate mean-field from beyond-mean- 
field effects, we follow the first strategy. Note that no isospin symmetry breaking effects 
due to Coulomb or charge-independence breaking of the nuclear part of the interaction 
is considered in the present study. It was shown to have no effect on the gap in the ^Sq 
channel [56]. Also, we do not consider the neutron-proton component of the force in this 
channel. 

2.2 Form of V^\f 

As already said, the bare A^A^ interaction has to be considered in the pairing channel at 
lowest order in irreducible vertices. However, the full complexity of any realistic NN force 
makes systematic HFB calculations in finite nuclei untractable from the computational 
point of view. We thus have to define a simplified bare force retaining the essential 
physics provided by the full NN interaction as regards pairing. 

A particular feature of the NN force in the ■'^S'o channel is the corresponding very large, 
negative scattering length. The empirical values for neutron-neutron, neutron-proton and 
proton-proton scattering length are a^'^° = — 18.5ib0.4 fm [57] (—18.7 ±0.6 fm in a recent 
experiment [58]), a^fo = -23.749 ± 0.008 fm [59] and a^fo = -7.8063 ± 0.0026 fm [60] 
respectively. This indicates that the NN interaction holds a virtual state in the vacuum 
at almost zero scattering energy in the ^Sq channel. In the vicinity of the virtual state, the 
scattering t-matrix can be written in a separable form, suggesting that the bare interaction 
itself is to a good approximation separable and non-local at low energy [61]. Thus, we 
start from the definition of the interaction Vg^° in the plane-wave basis 

^k,sS^= {r\kiSi) =e^^>^xi^ , (3) 

through 



( kiSik2S2 I V^gpO I ^353^454)0 = 2 ^ ^1^2 I ^sep° I 4^4)5 {SsiS3 ^3284-^8184 ^8283) 
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ihhlK^l'^lhh) = {k\V,\l'>\k'){2nf6{P-P') , 
where its center of mass part is approximated by: 

{k\K'f^°\k') = Xvik)v{k') . (5) 

While ki denotes the momentum of a particle in the laboratory frame, P = ki + kj 
and k = {ki — kj)/2 are the total and relative momenta of a pair, respectively. The states 
{| kiSi)} span the tensor product of momentum and spin (sj = it^) single-particle Hilbert 
spaces and are orthonormalized through* 

ih^.sJ'Pk.s,) = {27rf S{k, - kj) 6,^s^ . (6) 

Also, the subscripts a and s in Eq. 4 denote antisymmctrizcd and symmetrized matrix 
elements respectively. A matrix clement with no subscript is neither symmetrized nor 
antisymmetrized. The isospin quantum number has not been specified since the form of 
the matrix elements in the {T = l,Tz = ±1} channels is trivial. 

2.3 Connection with scattering phase-shifts 

The link between the bare force and the NN phase-shifts is obtained by treating the two- 
body problem in the center of mass frame. To make this link, it is convenient to introduce 
the energy-dependent scattering t-matrix. The Lippmann-Schwinger equation [61] defining 
it in the uncoupled channel takes the form: 



{k\t'%s)\k') = {k\v'^<^\k') + J 0^{k\v''^\k")Fi^^„is){k"\t'<^{s)\k') . 

(7) 

where P is the conserved total momentum of the pair and 



-ky 



s — m-\- le 



is the non-interacting two-particles propagator in free-space appropriate for outgoing 

boundary conditions. In Eq. 8, m is the nucleon bare mass and s is the energy of the 

interacting pair in its center of mass (total energy subtracted by P^/4m). It is worth 

noting that is independent of P and diagonal in k. 

* Because of the convention used to define plane waves, integrals in momentum space are characterized 
by / d^k I (27r)^. We also use the convention ti" = 1. 
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Because of the rank-one separable form chosen for the bare force, the Lippniann- 
Schwinger equation is exactly solvable and t''^'^ takes for any triplet (s, k, k') the separable 
form [61]: 



{k\t^°{s)\k') = \v{k)v{k')/ 



, dk" k"'^v^k") 

1 — A 



/ 

Jo 



(9) 



poo 

tan (5 ^° (k) = -Xrmrkv'^{k) / tt^ - 2Am V / dk" 

I Jo 



iQ 27r^ s — k"^/m + i€ 
While the full scattering problem is expressible in terms of the half on-shell t-matrix, 
the phase shifts carrying the information about the two-body wave-function at long dis- 
tances relate to the fully on-shell part of t [s = k'^/m = k'"^ /m). With the convention 
chosen to define the plane-wave basis, this link can be written explicitely under the form: 

{k\t{k^/m)\k') = -i^ V(2Z + l)e*'^'(^)sin5,(A:)Pi(cosA;iO , (10) 

I 

where 5i{k) denotes the phase shifts for a relative orbital angular-momentum I. We have 
not considered in Eq. 10 the coupling between different I channels as provided by the 
tensor force. Focusing on Z = and using the separable f-matrix, one gets: 

k"^v\k") ] 

^2 _ y^// 2 ' V / 

where V denotes the principal value of the integral. 

Since the separable form of the bare interaction in the ^Sq channel is physically moti- 
vated, one can hope to define an efficient and simple enough interaction by plugging the 
known phase-shifts into Eq. 11 to fix its parameters. Indeed, it has been shown that the 
^5*0 gap in nuclear matter is entirely determined by the NN scattering phase-shifts in the 
vacuum [56, 62]. Of course, a rank-one separable interaction cannot reproduce the phase 
shifts up to infinite energy in this channel since they change sign around 250AIeV [56, 63]. 
One could use a rank-two separable form to take care of this [63]. However, an overall 
reproduction of the phase shifts up to Eiab = 250 MeV <^ k = 1.73 fm~^ should be suffi- 
cient to describe pairing at relatively low density. Indeed, the kernel of the gap equation 
is strongly peaked at kp. 

The inverse scattering problem, which corresponds to the determination of a two- 
particle potential from the knowledge of the phase shifts at all energies, is exactly and 
uniquely solvable for rank-one separable potentials [61]. Thus, given the phase shifts, a 
unique solution exists for v{k). We do not proceed however through the resolution of 
the inverse scattering problem. Indeed, a simple analytic form of v{k) will eventually be 
necessary to perform HFB calculations of finite nuclei. We choose a simple form for v{k) 
and try to reproduce S^^° {k) as well as possible. We consider the Gaussian form: 



8 



v{k) = e-"'^' , (12) 

where a is the second parameter of the force. This choice will be shown later to be 
appropriate. 

2.4 ^5*0 pairing gap 

We have introduced a simple force to mimic the realistic NN interaction in the ^Sq channel. 
Starting from this interaction defined in the vacuum, one can go to the medium and 
compute the pairing gap through the BCS gap equation [65]'^. This scheme corresponds to 
the lowest order in the Goldstone-Brueckner-Bogolyubov perturbation theory [10, 11], or 
in the Galitskii-Gorkov Green function method [8, 66], and defines a meaningful mean-field 
picture. The simplest medium to consider at this stage is infinite nuclear matter. Indeed, 
its translational invariancc strongly simplifies the treatment and avoids the additional 
effects associated with the finitcncss of the nucleus. Of course, it is to some extent a toy 
model, even if neutron stars can be considered as being closely connected with it. 

In infinite matter, the favored Bogolyubov transformation correlates pairs of nucleons 
with zero total-momenta^ . The gap equation written in the plane- wave basis is of the 
usual BCS form and reads in the ^Sq channel as: 

Aw-A. = -^~^.-(giei«')^. (13) 

where Ek = [e^ — l^)'^ + is a quasi-particle energy; being the in-medium on-shell 
single-particle energy associated with the state 0^ and /j, the chemical potential. To be 
consistent, /x should be calculated iteratively by constraining the density of the system, 
while Cfc should be defined after regularizing the repulsive core of the bare force and by 
taking the influence of pairing correlations explicitly into account. This would typically 
require the use of the on-shell self-energy F^^ computed from the Galitskii T-matrix or the 
Brueckner ^-matrix defined in the superfluid phase [10, 11, 70]. A corresponding mean- 
field scheme is depicted diagrammatically on Fig. 1. Such a procedure is very involved, 
especially if using modern realistic NN forces in all {S, T) channels. It becomes prohibitive 
when dealing with finite nuclei. Thus, approximate schemes making use of the Brueckner 

Solving the gap equation using a separable force was done as early as 1964 [64]. In that work, qualita- 
tively similar results as those we derive in infinite matter were obtained using Yukawa type interactions. 
*Some exceptions exist however. See Refs. [67, 68, 69). 
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Figure 1: Definition of the normal Fg^ and abnormal Aj.^ parts of the self-energy at lowest 
order in irreducible vertices. Combining these definitions with the usual expressions of the 
BCS occupation numbers u|. and v^.., obtained from the compensation of the so-called 
dangerous diagrams at lowest order, provides the BCS gap equation [9, 10, 11]. The ir- 
reducible kernel entering the definition of T^^ at lowest order is the ^-matrix summing 
particle-particle ladders. The two-body propagator characterizing the intermediate states 
in the ladder is the product of one-body mean-field Green functions defined in the su- 
perfluid system. The irreducible kernel entering the definition of A^^ at lowest order is 
the bare NN force. The reason why the ^-matrix or the T-matrix cannot enter the gap 
equation originates from the necessity to avoid double-counting when compensating for the 
dangerous diagrams. It requires the exclusion of isolated particle-particle and hole-hole 
intermediate states in those vacuum-to-pair diagrams involving an abnormal contraction. 

G [71] or the Feynman-Galitskii T [66] matrices in the normal fluid are usually considered 
in infinite matter. 

When performing extensive mean-field calculations of finite nuclei, a phenomenological 
effective vertex such as the Gogny [20] or the Skyrme [41, 72] force is employed in the 
particle-hole channel to approximate one of the previous in-medium matrices. They usually 
incorporate additional phenomenology (and to some extent the effect of higher order terms) 
by fitting some experimental data at the mean-field level. As we are not looking for refined 
calculations at this stage, and because we want to use a clean theoretical quantity to 
adjust the separable interaction, we simply insert free single-particle energies = k'^ /2m 
into Eq. 13. We also take ji to be equal to kp/2m, where kp = (Svr^p)^/^ is the Fermi 
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momentum of one kind of nuclcons in the free gas at the density p. 

Considering these approximations and inserting the separable interaction into Eq. 13, 
the solution of the gap equation takes the form Ajt = Aow(/c), with the gap at zero 
momentum Aq satisfying the equation [56]: 

dfe \k^v\k) _ ^^^^ 

^""^ y/(A;2 /2m - /2m)2 + Ag (fc) ' 

After solving the Eq. 14, the gap at the Fermi surface is obtained through A^^^ = 

l^Qv{kF)- 

2.5 Fit 

We perform a combined fit of the force on the neutron-neutron scattering phase-shifts and 
on the pairing gap in infinite matter provided by the modern ylFlS NN interaction [73]. 
Among other features, ^F18 fits the proton-proton and neutron-proton phase shifts up 
to 350 MeV, as well as neutron-neutron low-energy parameters (scattering length and 
effective range) in the ^Sq channel. Thus, the neutron-neutron phase-shifts that we use 
beyond the validity of the effective range approximation k > 0.2 /m^^) are state of 
the art theoretical predictions. 

Thus, Eqs. 11 and 14 are solved for several set of parameters (A, a). The results are 
compared with those derived from ^yi8 with the same e^. The gap A^^ is plotted on the 
right panel of Fig. 2 for several values of the range ^/a; the intensity A being chosen in 
order to obtain the maximum at SMeV. Indeed, this is a solid prediction from all modern 
realistic forces [12, 56, 74, 75]. Note that the ^Sq gap calculated with free kinetic energies 
is very similar for all modern forces [12, 56, 74, 75]. This is due to the fact that they all 
reproduce the phase shifts very accurately. One can see that the gap strongly depends 
on a. As the pairing gap probes the interaction in a very sensitive way, requiring the 
precise reproduction of Aj^^^ derived from AF18 allows little latitude and determines the 
parameters of the separable force quite uniquely. This is a nice feature. For the best set 
of parameters (A = — 840 MeV. fm^, a = \/0.52 fm), A^^. is reproduced almost perfectly 
up to the gap closure. In particular, the bump is obtained at the right density and energy. 
This is a non trivial result in view of the very simple form of our bare force. The success 
of the procedure comes back to the justification of its separable form and to the overall 
reproduction of the phase shifts. 

The phase shifts calculated using the same sets of parameters are compared to those 
predicted by ^FIB on the left panel of Fig. 2. The simplicity of the force used seems to 
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Figure 2: Left panel: comparison between the neutron- neutron Sq scattering phase-shifts 
obtained from and from the rank-one separable force for several sets of parameters 

(A, a). Right panel: same comparison for the ^Sq pairing gap. In both cases, the range of 
the force is varied while the intensity is chosen accordingly to obtain a maximum gap at 
3 MeV. 

be more critical as regards a precise reproduction of the phase shifts. Also, this quantity 
hardly constraints the parameters in an obvious way. This justifies the complementary use 
of the pairing gap to fit the interaction. Interestingly enough, the gap is well reproduced at 
very low density when using our best set of parameters (thick curve on both panels), while 
the phase shifts are poorly reproduced below k = 0.4 fm~^. In particular, not reproducing 
the large scattering length a^'^° does not seem to be a major problem to obtain excellent 
gaps, as long as 5^^'^{k) is correctly treated beyond k = 0.4 /m~"^ as seen in the left panel 
of Fig. 2. It is known that concentrating on the very low energy part of the phase shifts (= 
effective range approximation) allows a good reproduction of the gap up to kp = 0.5 fm^^ 
but fails badly beyond that point [56]. As long as the quasi-bound state exists in the 
channel to motivate the separable form of the force, it seems non-essential to obtain a 
precise value of its eigenenergy (related to al^^, that is to the slope of S^^'^{k) at A; = 0). It 
seems more important to get an overall reproduction of 6^^°{k) which, through the on-shell 
t-matrix, relates to the wave- function of the virtual state at intermediate distances [61]. 
This result agrees with the conclusions of Ref. [74] and balances those obtained in the 
relativistic context [18]. 
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2.6 Analysis of V^\f 

As already mentioned, the rank-one separable form does not allow the description of the 
negative part of the ^Sq phase shifts at high energy. As seen on the left panel of Fig. 2, this 
translates into an overshoot of the phase shifts beyond k = 1.4 fm^^, but only into a slight 
over estimation of the gap at kp = 1.4 /m^^. Roughly speaking, the justification for not 
resolving explicitly the hard core of the realistic bare force is similar to the one providing 
the grounds for the V/q^ ^ interaction [76] or motivating the description of nuclear matter 
through effective field theories [77, 78]. Of course, we only take care of the ^5*0 channel 
here, we do not integrate out the high relative momentum components of the bare force 
explicitly and thus, we do not look for a high precision potential model. In Fig. 3, the 
matrix elements of our separable interaction in momentum space are compared with those 
of Ay 18 and those of the non-hcrmitian Vioy^k obtained from Ay 18 [79]. The cut-off used 
for Viowk is A = 2.1 fm^^ . Our y^pp" is very close to Vi^wk c'-iid quite different from ^41^18 
itself^. It clarifies the physical content of our separable interaction and characterizes our 
scheme as a low-energy effective theory of nuclear matter valid below ~ 3 Psat- Note that, 
while the matrix elements of modern realistic forces are scattered, their Viowk partners all 
look the same [79]. In the same way, the previous fitting procedure would lead to very 
similar separable interactions by starting from other modern realistic forces. 

As seen in Fig. 3, the range and the non locality of V^J^" arc of very good quality 
while its intensity at small k characterized by A = — 840 MeV.fm^ is slightly too low. 
It originates from the poor description of the phase-shifts below k = 0.4 /m~^, that is, 
from missing the virtual state at almost zero energy. This could not be avoided because 
of the simple Gaussian form used for v{k). Note that keeping the range fixed, a pole 
at zero energy would have been obtained in the t-matrix derived from F^g^" for A = 
— 950 MeF./m^. Last but not least, the ^Sq gaps obtained in infinite matter with our force 
and with Viowk are quite similar [80]. However, as pairing is exponentially sensitive to the 
force intensity, the slightly stronger diagonal matrix elements of Viowk up to A; = 1.5 fm~^ 
translate into a slight overshoot of the gap provided by AVIS or by our force [80]. 

The Gogny interaction can be considered as a benchmark concerning pairing properties 
in finite nuclei [6, 20, 42, 82, 83, 84, 85]. As regards the calculation of the ^Sq gap in infinite 
matter, it was shown to behave almost like a bare force, especially when using the DIS [81] 
parametrisation [38, 39]. To address more precisely whether the Gogny force mimics the 

^Vsgp° and Viowk differ from j4F18 by a constant shift in momentum space equivalent to a contact term 
in coordinate space. This contact term properly deals with the short range part of the NN interaction at 
low energy. 
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Figure 3: Matrix elements in momentum space of (full line), Viowk (dashed line), 

our separable force (dotted line) and the Gogny force with the DIS parametrization [81] 
(dashed-dotted line) in the ^Sq channel. Left panel: diagonal matrix elements. Right 
panel: non-diagonal matrix elements at k' = 0.009 /m~^. 

bare force in the channel, its diagonal and non-diagonal matrix elements [86] are also 
plotted on Fig. 3 for the DIS parametrization [81]. The Gogny force appears to be similar 
to the Viowk interaction and to our separable bare force. The underestimate at very low 
momenta is of no importance for pairing as discussed for Vg^° . The similarity with Viow k 
and Vg^° seems to explain why the Gogny force provides equivalent gaps to those obtained 
from the bare force. Looking more into detail, a slight overshoot of the gap provided by 
Viowk, and thus by AV18, is obtained with Gogny beyond kp = 0.8/m~^ [38, 80]. This 
can be related to its slightly too large diagonal matrix elements compared to those of our 
separable force beyond k = 0.75 fm~^. As densities beyond kp = 0.8 fm~^ dominate in 
finite nuclei, such an overshoot could be sufficient to explain why the bare force would 
not provide enough pairing in these systems and characterize the necessary beyond-mean- 
field effects. We will come back to this in section 4. Let us mention that we found 
the matrix elements of the Gogny force in the channel to be very similar to those 
of Viowk- Consequently, the Gogny force does not provide any artificial pairing in that 
partial wave. Such a good agreement in the D-wave was not expected since the Gogny 
force was, as any phenomenological force in finite nuclei, adjusted without taking great 
care of its partial- wave content. 

Due to the previous comparison, we have enough confidence in V^ep° to consider 
more detailed pairing properties. Beyond the ability of our force to reproduce A^^^, 
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Figure 4: obtained from AF18 and from the rank-one separable force for three different 
densities. 



it is worth analysing the momentum dependence of the gap at fixed kp. In Fig. 4, 
the Afc obtained from ^yi8 and from our force are compared for three different den- 
sities [87]. The agreement is excellent in the energy and density intervals of interest 
(A; < 1.5 fm~^, kp < /c|?* = 1.33 /m^^). The rank-one separable force is not designed to 
reproduce Ajt for k > 2 fm^^ where the NN phase-shifts become negative. However, this 
is not a significant problem to describe pairing at low energy and density. Note that the 
momentum dependence of the gap is quite insensitive to the realistic interaction used, at 
least at low energy [12, 74]. 

The good reproduction of shows the ability of the force to describe fine pairing 
properties. It confirms that our choice for v{k) is appropriate and that its range has been 
properly fitted. This comes back to the fact the gap Afc = AQv{k) is directly determined 
by the half on-shell t-matrix at s = (oc v{k)), that is by the vertex function of the virtual 
state in the channel [18]. 



3 Effective pairing interaction at the mean-field level 
3.1 Formalism 

From the previous discussion, the property of the Gogny force as being close to a reduction 
of the bare force for low-energy phenomena can be understood. On the contrary, in spite 
of their success as phenomenological pairing interactions, DDDI cannot be interpreted 
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as direct approximations of the bare force. One needs to understand them as effective 
vertices, which requires the derivation of an appropriate scheme. Also, our Vg^" bare 
force is still too complicated to be used in coordinate space HFB calculations. To deal 
with these two issues, we now recast the pairing problem in a slightly different manner. 

Let us start from the gap equation written at lowest-order in a given single-particle 
basis {aj}: 

^ij = ^{ij\'^\'>T^^) i^mn = ^ {i j IV \ mn) F^^g^{0) Aqr , (15) 
mn mnqr 

where k,... are convenient labels to denote that single-particle basis and V is the bare 
NN force. We define A and k as the energy-independent pairing field and anti-symmetric 
pairing tensor respectively (Kmn = ( ^ I I ^ ) > where | $ ) is the unperturbed quasi- 

particle HFB vacuum). For future convenience we have also introduced a two-body prop- 
agator through Kmn = T,qr ^mnqri^) \r- 

Let us now introduce a two-body amplitude TZ in the medium through the equation^: 

{ij\n{s)\kl) = {ij\V\kl) + J2 {ijms)\mn)FZ^,ris){qr\V\kl) , (16) 

mnqr 

where s is an external energy parameter and F^ a two-body propagator which will be 
specified later. Combining Eqs. 16 and 15, one can write: 

Y,{ijms)\kl)F,%{s)A,r = J2 iij\ns)\kl)F,l,{s){qr\V\st)F,t,{0)A, 

klqr klqrmnst 

(17) 

which, combined once more, allows recasting the gap equation under the form: 

A,,- = ^ ( i j I 7^(s) \mn) [F^,,,(0) - F^,^,(s)] A,, . (18) 

mnqr 

No approximation has been done to derive Eq. 18 from Eq. 15, and the derivation is 

valid whatever the two-body propagator F^ is. The solution of the pairing problem has 

to be the same by either using Eq. 15 or the set of coupled Eqs. 16 and 18. The previous 

rewriting of the gap equation is not restricted to the mean-field level. The propagators 

F^ and F^ may include off-shell nucleon propagation (fully or through a quasi-particle 

^We make use of discrete sums in Eq. 16. Of course, an appropriate integral is to be considered when 
dealing with continuous quantum numbers. The specification to the ^5*0 channel of the interaction is 
implicit in the present work even if Eq. 16 is valid in the general case. 
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approximation) at finite temperature [88, 90]. Eq. 18 still holds in that case. As we 
are interested in the mean-field treatment of the system at zero-temperature, we will only 
propagate the nucleons on-shell. Including dispersive effects in both the single-particle self- 
energy and the gap equation [70] constitutes a more advanced treatment of the many-body 
problem which is not conceivable nowadays for extensive calculations of finite nuclei. We 
have not tried to proceed to the same recast by starting from an irreducible vertex beyond 
the bare force. This could be done, in particular if using a quasi-particle approximation 
and/or considering a static limit for the higher-order terms in the interaction. 

The 7?.-matrix is specified as soon as the choice of the two-body propagator is 
made. Several choices are possible and selecting a particular one is simply a matter of 
convenience and of formal consistency. For several reasons discussed below, we require 
that the two-body effective vertex TZ treats particles and holes in a symmetric way and 
takes the superfluidity of the medium into account. More specifically, we define in 
such a way that TZ sums through Eq. 16 particle-particle and hole-hole ladders in the 
presence of pairing correlations. 

By considering the canonical basis associated with the actual Bogolyubov transforma- 
tion solution of the problem and the hypothesis that time reversal symmetry is conservedl' , 
is) and (s) read as [10, 11, 70]**: 



^ (19) 

^mnqris) = „ , cy ^rnq ^nr ^nm t 



where p„j = pmm = [1 ~ (^m ~ /^)/^m]/2 = Prn embodies the (diagonal) one-body density 
matrix and E„i = Em = a/ (cm — pY + ^mrri ^ quasi-particle energy. The indices (m, m) 
characterize the paired states in the canonical basis. The diagonal matrix elements of the 
single-particle mean-field e„j involved in E„i are defined using an appropriate scheme as 
already discussed in section 2.4. Because of the hypothesis wc have made concerning the 
Bogolyubov transformation, only the diagonal matrix elements A^j^j will be non-zero, 
selecting in Eqs. 18 and 19 the matrix elements of TZ and F^ involving paired two-body 
states of the same sort. 



"To be more specific, expressions 19 are rigorously valid only if the third part of the Bogolyubov 

transformation [89] is trivial. 

** Including off-shell propagation through a quasi-particle approximation would provide an additional 
factor ZjnZn, where Zi is the quEisi-particle strength [88]. 
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The first term entering the definition of -F„^„g,,. in Eq. 19 sums particie-particle (p-p) 
ladders as in the G-matrix except for the fact that BCS occupation numbers and quasi- 
particle energies appear because of the superfiuid nature of the system [10]. The same 
is true for the second term delaing with hole-hole (h-h) ladders. Two different signs 
are considered to sum the diagrams associated with h-h ladders in Tl. The "— " sign 
corresponds to what can be denoted as the T-matrix in the superfiuid system since Tl 
reduces in that case to the usual Galitskii T-matrix in the normal phase. This choice 
corresponds to summing the h-h ladders as they actually appear in the expansion of the 
ground-state energy when taking abnormal contractions into account in the theory [10, 11]. 
Then, we will denote the in-medium matrix associated with the sign as the P-matrix. 
Note that even if the definition of this P-matrix through Eq. 16 is fully valid, it does not 
correspond to the summation of hole-hole ladder diagrams as they appear in the expansion 
of the ground-state energy of the system. 

In the gap equation, it is appropriate to use T/P fully off-shell at s = 0^^. Eventually, 
Eq. 18 simplifies for the two considered reaction matrices into: 

m 

= 

m 

Again, different choices could have been made in Eq. 16 for the two-body propagator, 
amounting to using other two-body vertices than the T/P-matrices. For instance, choos- 
ing the plane- wave basis, reducing Eq. 16 to the center of mass frame of the interacting 
pair and using Eq. 7 for the two-body propagator, Eq. 16 would reduce into the Lipmann- 
Schwinger equation defining the t-matrix. Our choice can in fact be seen as an extension 
to the superfiuid medium of the regularizing procedures consisting of re-expressing the 
gap equation in terms of the scattering length [91, 92, 93] or of the t-matrix in the vac- 
uum [53, 74]. The latter choices are of particular interest in the low density regime where 
they allow the derivation of analytical formula for the gap in terms of known physical 
quantities. Eqs. 20 are not as useful in the low-density limit since the off-shell T/P(0) 
match a^^° only after the superfluity has disappeared. 

^^In Eq. 16, a half or fully on-shell matrix element of TZ corresponds to 2s — Ek + Ei. The choice 
s = corresponds to particular off-shell matrix-elements since it cannot be obtained by inserting any 
single-particle energies tk and ez into 2s = Ek + Ei, as long as pairing correlations are present. Indeed, 
the quasi-particle energies are always different from zero in such a case. 



\V{0)\mrh)2p„ 



Ar 



2Er, 
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However, taking T/D^O) is a formally optimal choice when studying nuclei at low 
energy. Such systems cover densities ranging from to psat, where the S'-wave like-particle 
superfluidity evolves from the weak to the intermediate BCS regimes, before coming back 
to the weak coupling. Thus, it makes sense to use a two-body scattering matrix taking the 
varying density and superfluidity into account. Also, considering the in-medium vertex 
at the threshold s = makes the effective gap equation as given by Eq. 18 to have its 
simplest possible analytical forms (Eqs. 20). The latter property is partly due to the fact 
the T/P-matrices treat particles and holes on the same footing which is reasonable when 
dealing with pairing correlations. For that reason, it should be preferred to the Brueckner 
G-matrix [62], even if calculated in the superfiuid system [10]. One could also have used 
the T-matrix in the normal phase. However, this vertex presents, like the G-matrix, a 
pole at the threshold {'^ s = fi foi the standard definition of its energy dependence) 
which signals the appearance of a Cooper bound state in the medium [62, 94, 95, 96]. As 
the aim of the gap equation is precisely to take care of the correlations associated with 
existing Cooper pairs, it should be combined with a regular vertex summing two-body 
correlations in the medium, except for those related to the possible bound-state. In fact, 
the Bogolyubov transformation aims to remove the pole from the scattering amplitude by 
defining a non-singular T-matrix at the threshold (s = in the definition used above), 
while treating pairing correlations explicitly through the gap equation [10, 11, 70, 97]. 

An immediate by-product of the previous derivation is to show explicitly that none 
of the previously discussed in-medium matrices can be used in the gap equation without 
re- writing the latter accordingly [12]. Aside from further aspects, the purpose of our 
recasting procedure is to incorporate the virtual high-energy transitions which appear in 
the original gap equation into the in-medium interaction. Whilst doing so, the important 
pair scattering around the Fermi surface are treated explicitly through Eq. 20. This is 
reasonable since the gap equation is almost linear in the high momentum regime while 
it is highly non linear around kp- Note that another recasting procedure was used in 
Refs. [12, 98] to get stable solutions of the gap equation when solved with realistic bare 
NN interactions. By summing virtual transitions above a sharp energy cut-off, an effective 
pairing interaction acting in a valence space was defined [12, 99]. It was anticipated 
and shown that the microscopic effective force was close to the off-shell T-matrix at the 
threshold [12, 100]. This result is not surprising in view of the previous discussion. 

The present scheme can now be translated into the definition of a microscopic effective 
pairing interaction to be used in the standard gap equation. By comparing Eqs. 15 and 20, 
we introduce two versions of such an effective vertex whose matrix elements in the canonical 
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basis read as: 



{ii\V^^^ \mm) = {ii\v'^^°{0)\mfh) 2 , 

(21) 

{ii\V^^^ \ mfh) = {ii\T ^°{{))\mfh) 2 {1 — pm) pm ■ 

Each of these two versions includes a smooth cutoff as well as in-medium correlations. 
The cut-offs have appeared naturally when recasting the gap equation and will adapt self- 
consistently to the amount of pairing in the system. They do not have to be additionally 
chosen or optimized. They are measured with respect to p and not with respect to the 
bottom of the mean-field potential and thus, evolve with density. Note that using T has 
led to the appearance of a symmetric cut-off on both sides of the Fermi energy while P 
comes together with a single cut-off above the Fermi energy. Again, one should stress 
that the effective interactions defined by Eq. 21 result from a rearrangement of the gap 
equation and that no approximation has been done in the mean-field treatment of pairing 
and that no effect beyond that level of approximation has been included at this point. 

The two versions of V^^^ have been studied. In the following, wc only present the 
results obtained by using the D-matrix. The reason for that is because the P-matrix is 
slightly less sensitive to the self-energy effects when those are included in the calculation 
as will be discussed in section 3.2.4. However, any important conclusion drawn in the 
following for the P-matrix is valid for the T-matrix. 



3.2 Infinite matter 

3.2.1 Calculation of the P-matrix 

Let us use the previous scheme to attack the infinite matter problem. The plane wave 
basis corresponds to the canonical basis associated with the Bogolyubov transformation 
solution of the problem. Using transparent notations, the two-body propagator involved 
in the equation of the V matrix reads as: 

FV^is) - _lZfMZff2zL (22) 

Performing the averaging over the angle (P, k) 

Fp,k{s) = ^ r r Fp^fe,eose(^)sin^d0# = /' dxF^,^Js) , (23) 

47r Jj^^oJg^Q J^^Q 
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and considering the separable form of Vg 
T>^^° in the center of mass yields: 




can be integrated explicitly. Expressing 



{k\v'^°{kF,P,s)\k') 



\v{k)v{k') 



(24) 



k"^v^{k")Fp^^„ {s) 



= Xv{k) hikF,P,s) v{k') . 



(25) 



Starting from our separable force, Eq. 24 shows that the in-medium vertex T> o is 
also separable in the three variables {k,P,k'). One has to check that using P^'^" in the 
gap equation as given by Eq. 20 provides the same gap as obtained in section 2.4. Since 
only pairs with a zero total-momentum occur in nuclear matter and since the integration 
of P for P = does not require the angle averaging procedure, this reproduction must be 
exact. We have checked that it is so, both for the gap at the Fermi energy Ai^^ and for 
Afc. In the following, we want to study and parametrize the function h characterizing the 
dependence of P^"^" on kp and P. This will be necessary to use V^-^-^ in HFB calculations 
of finite nuclei. 

3.2.2 Density dependence of V^^o (q) 

We now study the density dependence of P^'^" (0) in infinite matter. It enters the effective 
pairing force as a factor defined as C{kp) = h^kp, 0,0). 

This function, derived using V^ep°5 Eqs. 24-25 and free single-particle energies, is shown 
in the left panel of Fig. 5. The many-body effects in P^"^" are such that the magnitude 
of the in-medium interaction decreases with increasing density and saturates for kp > 
lfm~^. Within the usual framework of the local density approximation (LDA), this 
property translates into a pairing force which is slightly enhanced at the surface of the 
nucleus kp ~ 0.8 — 1.1 fm,^^) as compared to the center (44> A;|?* ~ 1.33 /rn^^). 

In infinite matter, the density dependence of the force has of course no link whatsoever 
with a surface effect. It is a pure density effect generated by the in-medium coupling of 
mean-field particles. It is only through the language of the LDA that one would talk 
about a position- dependent interaction in the nucleus. The present study clarifies this 
character of the effective interaction at the mean-field level. Additional effects, coming for 
instance from the induced interaction generated by the exchange of surface vibrations in 
finite nuclei, have to be considered on top of the mean-field approximation [55]. 

The diagonal matrix elements of V-^^-^ at kp are compared with those of Viowk and of 
our separable force in Fig. 6. Generally speaking, the pairing gap A;^.^ is not determined 
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Figure 5: Left panel: derived density dependence of the effective pairing force in infinite 
matter (full line). Two different fits are also given. Right panel: same for the zero-range 
approximation of the effective pairing vertex. 

by the diagonal matrix elements of the force at kp only [74]. Otherwise, the effective 
interaction, our separable bare force and AV18 could hardly give the same gaps. Obtaining 
the same gaps is possible because the different nature of their diagonal matrix-elements 
is compensated by different off-diagonal characters. For instance, the reduced influence 
of the off-diagonal processes in the case of the effective force embodied by the cut-off 2 pk 
is accompanied by an enhancement of its diagonal matrix elements through the density 
dependence. This shows why the in-medium effects resummed in the effective pairing force 
are correlated with the the cutoff emerging through the recast of the gap equation. Note 
that only in the case of such an effective interaction resumming off-diagonal processes (up 
to high energy if starting from a realistic bare A^A^ force) the gap is indeed very much 
determined by the matrix elements at kp. This would be particularly true if using T 
instead of D, since then the off-diagonal matrix elements are cut on both sides of the 
Fermi energy as shown by Eq. 20. Only this property authorizes the use of weak coupling 
formulas [52]. 

Later on, we will need C{kp) as an analytical function of kp- For that purpose, we 
fit C{kp). One should stress that such a fitting procedure does not correspond to the 
use of free parameters in the definition of the pairing interaction. We simply reproduce 
/i(A;i?, 0, 0), derived with no freedom from Vg^p , itself fitted on scattering properties and 
on the ^5*0 pairing gap provided by ^FIS. As the gap is exponentially sensitive to the 
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Figure 6: Diagonal matrix elements of different interactions at kp in the Sq channel as 
a function of kp in infinite matter: Vigyjk (full line), the separable bare force (dashed- 
dotted line), the finite-range effective pairing interaction (dashed line) and its zero-range 
approximation (dotted line). 

strength of the interaction, a fine fit of C{kF) is required to reproduce A^^^ with high 
accuracy. Two different fits are compared with the exact function C{kF) on the left panel 
of Fig. 5. The gap re-calculated using the fit defined through powers of Inkp is compared 
with the gap derived from the bare force on Fig. 7. They are nearly identical. 

The two fits displayed on Fig. 5 have a quite different analytical character. Using 
enough terms, one can actually propose a large set of precise fits making use of very 
different power series. For a given number of terms, some expansions are more precise 
than others. For instance, an expansion in powers of kp requires a large number of terms 
whereas expansions in powers of \/k,p or kp^^^ . . . converge much faster. Whatever the 
chosen expansion, the important feature is to reproduce precisely the behavior of C{kp) 
derived from V^^'^. 

The fit making use of powers of In kp provides an excellent approximation of C{kp) 
with only three terms and presents the particular feature that adding other powers of In fc^? 
does not improve the fit significantly. As this was not the case for any other expansion 
we have tried, we consider this function as an "exact" analytical form of C{kp). This 
function rises in the limit kp 0. This agrees with the fact that P^'^"(0) tends to the 
off-shell t^'^"-matrix at s = 0, which almost diverges because of the virtual state in the 
vacuum at approximately zero scattering energy. This is seen in Fig. 6. As no real bound 
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Figure 7: ^Sq pairing gaps obtained from the bare force (full line), the effective force 
(dashed line) and its zero-range approximation with two different fits of its density depen- 
dence (dotted and dotted-dashed lines). 

state exists in this channel, the pairing collapses accordingly in the kp ^ limit as seen 
in Fig. 7. 

3.2.3 Zero-range approximation 

An appealing property of the scheme developed here is its suitability for studying a zero- 
range approximation of the effective pairing interaction. Since the high-energy transitions 
appearing in the gap equation have been resummed into the in- medium vertex 7?^^°, no 
divergence is expected when taking its zero-range limit. This is embodied by the presence 
of the cut-off 2 pk in the gap equation. 

Expanding {k\'D^^°{kF, P,0) \k' ) in the range and in the non-locality leads to an 
effective pairing interaction reading in momentum space as: 

( k I V'ikF) I ' ) « 2 A C'^ikp) Pk' ■ (26) 

One is simply left with a density-dependent force, constant in k and varying as pk', 
where k' is the momentum of the intermediate state in the gap equation. The local 
approximation consisting of taking V^'^ independent of P is shown to be appropriate in 
section 3.2.5. The gap equation in infinite matter associated with the force defined through 
Eq. 26 provides a constant solution A = AfT . 



24 



The structure of the zero-range approximation of the effective pairing force we obtain 
here provides a formal justification for the use of DDDI complemented with a smooth 
cut-ofF above the Fermi energy [102]. Considering the T-matrix instead of the P-matrix 
would allow a similar justification for those DDDI which have been used with a symmetric 
cut-ofF on each side of the Fermi energy [22]. 

We now extract the function C^*'(A;f). This is done by matching A^^ with A^^ ob- 
tained from AV18 as a function of kp- The function C^'^^kp) is plotted on the right panel of 
Fig. 5. As expected, no difference is seen at very low density between C{kp) and C^^^kp), 
whereas the effect of the finite range becomes more and more important with increas- 
ing density. While C{kF) was almost fiat and non-zero when approaching the saturation 
density of nuclear matter, the zero-range in-medium pairing interaction goes to zero. In 
other words, the present calculation predicts (through the LDA) that, when approximated 
by a zero-range-like vertex, the effective pairing interaction has to be renormalized by a 
density-dependent intensity whose surface character is more pronounced than in the case 
of the finite-range original version. Quantitatively speaking, while Cikp) is multiplied by 
a factor 1.4 between k^-* = 1.33 /m~^ and kp = 0.8 fm~^ {■^ psat/^), C^'^ikp) is multi- 
plied by a factor 3.6. The zero-range vertex has a behavior between surface and volume. 
This confirms the results obtained from refined phenomenological studies performed with 
these kinds of vertices [23, 30, 31, 32]. 

The different behaviors of the finite-range and zero-range forces can be understood in 
the following way. The density dependence saturates when the size of the Fermi sea is 
of the order of the inverse of the interaction range. Beyond that point the range governs 
the coupling inside the Fermi sea. This is visible on the left panel of Fig. 5 where the 
finite-range interaction saturates beyond kp > l/\/a ~ 1.4 /m~^). For the zero-range 
vertex the same only happens when kp ^ l/^fo. = oo where the interaction goes to zero 
in order to compensate for its artificial constant coupling in A;-space. 

It becomes clear in the present context that the surface enhancement of a zero-range 
pairing vertex takes care to some extent of the finite range of the nucleon-nucleon in- 
teraction. The effect of the range, although noticeable, seems to be re-normalizable at 
the mean-field level. This is seen on Fig. 6 where the diagonal matrix elements of the 
two interactions at the Fermi surface are compared. They are nearly identical. The 
possible renormalization of the range conveys that nuclear matter presents a so-called 
weak-coupling BCS regime over the density range of interest characterized by the size of 
the Cooper pairs being much larger than the inter-particle distance and the range of the 
force. Such a result provides the grounds to the local density functional theory [14, 15] in 
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the context of superfluid nuclear matter [44, 45], even if the context is shghtly different 
in that case since all beyond-mean-field effects are included in the in-medium coupling 
constants. 

However, the fact that the range and the non-locality of the force do not need to 
be resolved has to be confirmed in finite nuclei. Indeed, a zero-range vertex carries less 
information than the finite-range force does; this being embodied in the present case by 
the different non-diagonal matrix elements of our finite- and zero-range forces. While the 
zero-range version is able to reproduce the gap at the Fermi level as a function of the 
density, it predicts constant gaps as a function of momentum at a given density. Such an 
approximation is doubtful, especially around the Fermi energy where the gap is rapidly 
varying [12], as seen in Fig. 4. Treating the range and the non locality of the force affects 
the state-dependent pairing gaps, the particle-width of deep-hole states in finite nuclei 
and translates to some extent into the spatial character of the pairing field [42]. It may 
also be of importance to describe excited states in nuclei. Resolving this issue requires an 
excursion far from the valley of /3 stability [42] . 

Using C^^{kF) within the LDA requires the identification of its analytical dependence 
on kp- Unlike the function Cikp), derived from the bare force through the calculation of 
the P-matrix, fitting C^^{kF) amounts to fixing the free parameters entering the definition 
of the phenomenological zero-range vertex. Two examples are plotted on the right panel 
of Fig. 5. The overall multiplication by A in Eq. 26 does not correspond to adding another 
free parameter. No parameter is needed to specify the cut-off. 

Fig. 7 displays the gap re-calculated by inserting the two used fits into Eq. 26. The 
three-parameter fit allows a satisfactory reproduction of the gaps derived from the finite- 
range interaction. The fit employing two parameters does not provide sufficiently precise 
results. This last parametrization overshoots the gap around its maximum by half an MeV 
and undershoots it beyond kp = 1.25 fm~^. 

The number of parameters usually used in connection with zero-range forces is two for 
a pure- volume pairing and three for a density-dependent pairing [21, 22, 38, 39, 101]. The 
reproduction of the gap in nuclear matter [38, 39] as well as recent calculations in finite 
nuclei [101, 23] favor a density dependence. In such studies, the exponent of the density 
has sometimes been varied, in which case it could be considered as an additional free pa- 
rameter. Even if the presently derived zero-range force also makes use of three parameters, 
it corresponds to a cleaner parametrization since the cut-off appeared naturally and the 
adjustment of the force is simply a matter of reproducing a fixed, derived function. 

Following the path of a range expansion, one could include correction terms in k^,k'^,k^* 
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k'"^ . . . and obtain the corresponding density dependence by using the previous method. 
Finally, it is worth noting that the finite-range effective interaction makes use of only two 
parameters for an even better quality of results. 



3.2.4 On-shell self-energy and isospin dependence of T>^^° (0) 



The logarithmic density dependence of the effective interaction has been derived from the 
integration of the P-matrix to all orders in the bare interaction. This integration includes 
both hole-hole and particle-particle ladders in the superfiuid phase and treats the energy 
denominator of the two-body propagator explicitly. However, while everything has been 
consistently compared so far using free single-particle energies, the parametrizations of the 
functions Cikp) and C''^^{kp) obtained without including self-energy effects are of course 
approximate. As mean-field calculations of finite nuclei will eventually be performed, the 
function C{kp) should be derived accordingly. Consequently, we recalculate the effective 
interaction including on-shell self-energy effects through a density-dependent effective- 
mass approximation for the single-particle energies: 



where 4 = 37rV/2 = [{k^f + {l^Ff]/2 and /? = - Pp)/p = [{k^f - (^^)3]/24 are 



the Fermi momentum and the matter asymmetry, respectively. Quantities relating to a 
particular isospin value present an index q either specified as n and p or by 1/2 and —1/2 
when dealing with neutron and proton, respectively. 

We use the effective mass as given by a standard paramctrization of the Skyrme 
force [103]. Such an effective mass mimics an average of the Brueckner-Hartree-Fock 
(BHF) fc-mass: 



over the Fermi sea. It is independent of fc, smaller than one at all densities and of the 
order of 0.7 at saturation density of symmetric nuclear matter [104]. The Skyrme effective- 
masses in symmetric and neutron matter are plotted on the left panel of Fig. 8 for the SLy4 
paramctrization [103]. As seen, self-energy effects are larger in symmetric matter than in 
neutron matter at a given density because of the stronger neutron-proton interaction. 

As the gap equation is particularly sensitive to the density of states at the Fermi 
surface, the use of an averaged fc-mass which does not reflect the bump at the Fermi 
surface of the actual BHF fc-mass is questionable. When dealing with pairing, one should 



2m*{kF,P) 



(27) 




(28) 
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Figure 8: Left panel: neutron effective mass in neutron (dashed line) and symmetric (full 
line) matter for the SLy4 Skyrme force. Right panel: ^Sq gaps calculated with three choices 
of single-particle energies: free kinetic energies (full line), effective mass approximation 
in neutron matter (dashed line) and effective mass approximation in symmetric matter 
(dotted-dashed line). 

maybe consider the BHF A;-mass at kp instead [105], or the complete BHF single-particle 
energies. However, two arguments are in favor of the averaged effective masses. First, 
the bump of the A;-mass at kp is smoothed out to some extent by the presence of pairing 
correlations [106]. Second, the main purpose of the present section is to include self-energy 
effects in V^^° which should not be as sensitive as the gap to the density of states at kp- 

The corresponding gaps for symmetric and neutron matter are compared on the right 
panel of Fig. 8 to the one obtained using free kinetic energies. The gaps are plotted as a 
function of the Fermi momentum of the species (neutrons or protons) concerned by the 
pairing. As expected, the inclusion of on-shell self-energies reduces the gap. However, 
the effective mass obtained from the SLy4 parametrization is such that the reduction of 
the gap in neutron matter is too strong compared to the one obtained using realistic bare 
interactions and BHF single-particle energies [12, 107]. On the other hand, the reduction is 
slightly too small in symmetric matter, especially at low density [12, 107]. This corresponds 
to isoscalar and isovector parts of the effective mass which are respectively too large and 
too small. In any case, the global effect is present at the densities of interest. This should 
be sufficient to discuss the effects of the self-energy on the effective pairing interaction. 

The modified functions C{kp) and C^^{kp) have been re-calculated for neutron and 
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symmetric matter. They are nearly identical to those derived using free kinetic energies, 
at least up to saturation density. Therefore, one can safely state that the vertex only 
depends on the density of the interacting nucleons Pq = kp/Sir^ = p[l + (— 1)^/^^'' /3]/2 
and is not influenced by the surrounding nucleons of the other species. This result offers 
a microscopic answer to the issue dealing with the isovector density dependence of the 
effective pairing force at the mean-field level. A dependence on the total density as used 
so far is not justified from the microscopic point of view. This could of course change 
beyond the mean-field [53]. 

At the present stage, the best logarithmic fit of the density dependence of the finite 
range effective pairing interaction is given when using three terms by: 

Cikp) « 0.978444 - 0.682204 In A;^ + 0.761575 (In kpf , (29) 
while in the case of its zero-range approximation, we find: 

C"^''(fcF) f« 0.420637 - 1.012900 In feiT + 0.708922 (In /ciT')^ . (30) 

Using these fits with the Skyrme effective mass approximation, the gaps obtained in 
neutron and symmetric matter reproduce well those shown in Fig. 8. Finally, note that the 
energy functional of normal and abnormal densities obtained with such density-dependent 
vertices in the particle-particle {T = l,Tz = ±1) channels is isospin symmetric. 

3.2.5 Total-momentum dependence of 'D^'^" (0) 

We now study the dependence of the P-matrix on the total momentum P of the interacting 
nucleons. As nucleons with non-zero total momenta can be paired in finite systems, it is 
important to get insights into the corresponding component of the effective force. In order 
to obtain the P-matrix for P 7^ 0, a numerical averaging over the angle between the 
relative and total momenta of the pair in the intermediate states has been performed in 
Eq. 23. 

The function h{kp,P,0) is plotted on Fig. 9 as a function of P, up to P = 2.5 fm~^, 
for densities ranging from kp = 0.2 fm^^ to kp = l.Afm^^. The value P = 2.5 fm~^ 
corresponds to the maximum P a pair of nucleons inside the Fermi sea can have at kp = 
1.25 fm^^ . The P dependence of T> is significant at extremely low density, whereas it 
becomes less pronounced with increasing kp. Around saturation density, the interaction 
is almost independent of P. Also, the interaction is strongly modified at low P when 
increasing the density while it is less sensitive to the medium for a pair having a large 
total momentum. 
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Figure 9: Total momentum dependence of the P^'^°-matrix for several densities (full lines). 
A comparison is done with the function C{kp) e""^^^/^ (dashed lines). 

For a reason which will be clarified later, we need to test the hypothesis that h{kF, P, 0) 
is separable in (kp, P) and that the total momentum dependence can be parametrized by 
f{P) = e^"^^^/^, a being the same range as before. As a matter of comparison, the 
product C{kp) e""^^^/^ and the function h(kp, P, 0) are plotted on Fig. 9. The agreement 
is not satisfying at extremely low density where h decreases much faster than C{kp) f{P) 
as a function of P. While the parametrization works well for typical surface densities, it 
decreases too quickly as a function of P around saturation density. Of course, a pair of 
nucleons formed in finite nuclei has non-zero components for all values of P. In addition, 
the intensity of the interaction correlating the pair for each of these components must be 
seen as an average over all densities. As a result, we expect to have a good description of 
the combined density and total-momentum dependences of the effective pairing force by 
parametrizing h{kp,P,0) through C{kp) f{P). 

3.3 Effective pairing interaction in coordinate space 

To obtain the previous results in infinite matter, we started naturally from the bare force 
expressed in the plane-wave basis and derived the effective pairing interaction accordingly. 
To perform calculations in finite nuclei, and especially if solving the problem in coordinate 
space, it is necessary to have the expression of the interaction as a function of the inter- 
acting nucleon positions. As already mentioned, the force presently used is finite ranged 
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and non local. Consequently, it depends on the four position vectors {^i, r2, rj^, of the 
incoming and outgoing interacting nucleons. Also, the effective nature of the force has 
been characterized by its density dependence C{kp). As the spin-isospin part of the force 
is trivial, let us consider the spatial part of the two-body interaction under the form: 

(nr2|P;^o(0)|f3f4) = j^^^ I drC{p,{r)) e'^'- l^^-^^^l'/^"' , (31) 

where 

C{pq{r)) = 1.179341 + 0.345992 Inpg(f) + 0.084619 [lnpq(f )]^ , (32) 

is the form factor resulting from the local-density approximation of the function C{kp) 
derived in infinite matter and discussed in section 3.2.2. The fit expressed through Eq. 29 
and the relationship kp = (Stt^p^)^/^ has been used to write Eq. 32. 

The form given by Eq. 31 is reasonable since, if one forgets about the density depen- 
dence (or takes it to be constant in space), the following identity holds: 

{nf,\vl''^{0)\rsu) = (^,^3)3/2,9 • (33) 

Hence the interaction simply generalizes the Gaussian form to the non-local case. The 
matrix elements of the interaction defined by Eq. 31 calculated in the plane- wave basis 
take the form: 

{hhl (0) Ihh) = {2Trf \C{kl) e~"' l'^' + ^ + 5(P - P ') . (34) 

This is precisely the P^'^" -matrix derived in infinite matter from our separable bare 
force and studied in previous sections. To reach such a form, one single approximation 
dealing with the P dependence of X'g'^" was performed in section 3.2. Thus, Eq. 34 clarifies 
why the approximate dependence on P was needed in order to write the effective pairing 
interaction in coordinate space. Note that the parameters of the force have been fixed in 
infinite matter and that no room is left for any adjustment in finite systems. 

It is worth characterizing the physical content of the LDA used for the function C{kp) 
in Eq. 31. The particle-particle and hole-hole ladders (including all high energy processes) 
associated with two-body scattering in the medium are resummed in the effective interac- 
tion by considering the nucleus as slices of homogeneous, infinite nuclear matter. However, 
the scattering of pairs around the Fermi surface responsible for the formation of the bound 
states in the medium are treated explicitly in the finite system through the resolution of 
the HFB equations. Also, the quantum shell effects which strongly influence the latter 
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pair scatterings arc taken into account explicitly. All along the way, the finite range and 
the non-locality of the interaction is fully considered. We expect the local approximation 
for the correlations associated with two-body scattering to satisfactory and the inclusion 
of gradient terms^ to improve on it to be rather unnecessary [108]. The fact that the 
LDA is performed on the effective vertex itself makes the present approximation different 
from a so-called semi-classical treatment [109] or from the local density functional the- 
ory [14, 15, 44, 45]. In particular, it has allowed us to derive fine details of the interaction, 
such as its non-analytical low-density or isovector character, which are otherwize difficult 
to identify by looking directly at total-energy differences because of our presently limited 
experimental knowledge [110]. 

One of the most important properties of the newly defined pairing interaction can be 
identified by calculating its antisymmetrized matrix elements in any given single-particle 
basis. To be explicit, we define a basis {fnpCg} where {n,p,(,q) denote the principal 
quantum number, the parity, the z-signature and the isospin of the state, respectively. 
This is typical of the HF or canonical basis of a (triaxial, rotating, odd) deformed nucleus. 
Specifying (pnpCqs ^ the component of cp having a good projection s of the spin on the 
z-axis (deformation axis in the intrinsic frame), the antisymmetrized matrix elements of 
the effective pairing force as defined through Eqs.21 and 31 read as: 



iV,'")^kJl = A(p,+p;)^y dfC(p,(0) ^*^,^^^,,(rl^,\^^^^^^^ (35) 

ss' 

X {'fujp^Qgsi'f) VnipiCigs'ir) - (pnjp^Qqs'ir) (PnipiCiqs{r)} , 

where (pnpCgs is defined as: 

^npCgsir) = j^^^ J df' e-\^^-^^'\y'-\^,^,s{r') . (36) 

Restricting the pairing to states of opposite signatures in the Bogolyubov transforma- 
tion provides the matrix element 35 with an additional factor (5-^.^^. ^-^jCr Writing the 
matrix elements in the HF basis requires the use of the cut-off {pj + pi) instead of 2 in 
the canonical basis. This is a natural, but not fully rigorous extension as can be realized 
by going back to the recasting procedure proposed in section 3.1. One should preferably 
work in the canonical basis. 

Solving the HFB problem in coordinate space is a nice feature since it allows a natural 

treatment of all kinds of deformations and is well suited to describe exotic systems for 

•'■■'■This is of course to be differentiated from gradient terms simulating beyond- mean-field effects like the 
exchange of surface modes. 
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which asymptotic properties of individual wave- functions and densities must be considered 
carefully. However, the naive use of a finite-range interaction in this context is numerically 
prohibitive [5, 46, 111]. Indeed, solving the BCS gap equation is too costly for practical 
applications, while the HFB problem takes the form of a coupled set of integro-differential 
equations, also untractable for systematic studies. This is the main reason why zero-range 
forces have been extensively used so far [21, 54]. The matrix elements of our finite-range, 
non-local effective pairing-interaction look very similar to those obtained from a zero- 
range force [112]. The only additional cost is to replace the single-particle wave-functions 
tPnpCqs by their convoluted counterpart as defined through Eq. 36. This property makes 
the corresponding HFB problem in (3-dimensions) coordinate space tractable through the 
two-basis method [47, 112, 113], and almost equivalent computationally to the use of a 
local, zero-range force. In fact, only trivial modifications of HFB codes using this method 
are required. Adding a routine to convolute the HF wave-functions, using the proposed 
density-dependence and inserting the factor (pj+pi) when calculating the matrix elements 
of the pairing field in the HF basis is necessary. Of course, the pairing field provided by 
our force is non local. This feature prevents from solving the HFB problem through a 
direct diagonalization of the HFB matrix in coordinate space [46, 114]. 

It is also worth noting that the derivative of the pairing energy with respect to the 
density matrix {pji} when minimizing the total energy to obtain the HFB equations by 
must not be considered with the present force. Indeed, our effective pairing interaction 
is equivalent to the bare force in the pairing channel and its dependence on {pji} simply 
arises through a recasting of the gap equation originally written in terms of the bare 
force. Of course, this statement is no longer valid when going to higher orders, when 
renormalizing the effect of the three-body force [29] or when using the philosophy of the 
density functional theory. 

As an alternative to our finite-range and non-local microscopic effective vertex, one 
can use its zero-range approximation studied in section 3.2.3. When performing the zero- 
range limit appropriately, the convolution used in Eq. 36 becomes the identity operation. 
Consequently, the matrix elements of the zero-range vertex take the form given by Eq. 35 
with (pnipiC,iqs replaced by ^nmCiqs and C {pq{r)) replaced by the LDA of the function 
C"'{kp) parametrized accordingly. 

Last but not least, it is essential to note that the separable bare force defined through 
Eq. 5 docs not lead to simple calculations in coordinate space. Its finite range and partic- 
ular form of non-locality make it numerically Tintractablc in such a context. This is only 
by going to the in- medium vertex that one obtains an interaction as given by Eq. 31. 
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4 Discussion 



While very low-densities are not essential for stable nuclei, they become increasingly im- 
portant when going to exotic systems. Indeed, nuclei close to the neutron drip-line develop 
extended low-density halos or skins where pairing correlations play a crucial role. This is 
also true in neutron stars' crust. Thus, the behavior of the pairing force at low density 
is of great interest for the study of these systems. It is clear that usual DDDI miss the 
very low-density part of C^^(A;f) discussed in section 3.2.3. In Ref. [31], strongly increas- 
ing interactions at low density were simulated through small exponents 7 in Eq. 2. Such 
pairing interactions were disregarded, notably because of the unrealistic reduction of the 
two-neutron separation energy across the magic number N = 82. However, these interac- 
tions were used together with a phenomenological cut-off at a maximal fixed energy in the 
single-particle spectrum. In the present case, the cut-off 2 pm derived in connection with 
the density dependence will weight the low-density content of the force in a very different 
way. 

Also, the isovector character of the pairing force should manifest itself when studying 
drip-lines systems. The dependence of our force on the density of the interacting nucleons 
will provide a weaker pairing when going toward the drip-lines than standard DDDI de- 
pending on the total density and adjusted around the valley of stability. One can expect 
these properties of the force to significantly influence matter densities, pair densities, in- 
dividual excitation spectra, low-energy vibrational collective modes, rotational properties 
and the odd-even mass differences in exotic nuclei as well as the position of the neutron 
drip-line. It was shown that the pairing gap is extremely sensitive to the details of the 
force when dealing with neutron rich nuclei [31, 115]. Also, while the position of the proton 
drip-line appears to be quite robust, the position of the neutron drip-line was shown to 
be shifted by up to 25 mass units depending on whether the so-called volume or surface 
pairing force was used [115]. This is one of the achievements of the present work to propose 
a pairing interaction whose ab-initio derivation should make its extrapolated use to exotic 
systems reliable and subject to less uncertainty. 

At the mean-field level, this ab-initio character deals with the bare force. This is 
the first essential piece since, unlike for normal superconductors in condensed matter, 
the bare NN force provides pairing between the constituents of the nucleus. However, 
the question arises of its contribution compared to the pairing generated by collective 
effects. Recently, the BCS gap provided by the realistic ^F18 NN force was shown to 
account for only half of the experimental odd-even mass staggering in ^^'^Sn [55]. This was 
interpreted as a necessity to go beyond the mean-field to introduce the off-shell nucleon 
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propagation associated with the particle-vibration coupUng and the induced interaction 
generated by the exchange of the same surface vibrations between time-reversed states. 
Including those processes, the experimental odd-even mass staggering was reproduced 
in ^^°5n [55]. This exchange of surface vibrations should provide the effective pairing 
force with an additional surface-peaked character and some isospin dependence due to the 
appearance of new vibrational modes toward the neutron drip-line. On the other hand, 
the same category of diagrams is known to decrease the gap in the bulk [2, 49, 50, 51, 52]; 
even if the corresponding dressing of the vertex could change significantly when going from 
neutron matter to symmetric matter [53] . The latter would add up another dependence of 
the interaction on the nuclear asymmetry. Thus, the situation concerning the net influence 
of beyond-mean-field effects on pairing in finite nuclei is unclear. Let us repeat that, as 
already discussed in section analysis, calculations performed with the Gogny interaction 
enforce the idea that medium renormalizations on top of the bare force in the T = 1 
channel should not be that large when going to finite systems. 

However, one has to consider to more elements before concluding. First, the effect of 
three-body force on pairing should be treated if one works with microscopic interactions 
as proposed here. The three-body force has been shown to decrease the gap in infinite 
matter non-negligibly for kp > 0.8 fm~^ [107]. Including the three-body force in the 
pairing channel could also be necessary to reproduce delicate phenomena such as the 
odd-even staggering and the kinks of differential charge radii [29]. Second, one has to 
include the Coulomb force in the proton-proton pairing channel. Prom the comparison 
between proton-proton and neutron-neutron scattering phase-shifts in the ^Sq channel [73], 
one would not expect a strong anti-pairing effect from the Coulomb interaction. This is 
particularly true since the very low energy regime where the Coulomb force is active does 
not seem to matter so much when dealing with pairing as discussed in section 2.5. However, 
it was shown that the self-consistency of the HFB calculations makes the Coulomb anti- 
pairing effect quite significant in finite nuclei [116]. 

Thus, pairing correlations and details of the effective pairing interaction are far from 
being understood in nuclei. Systematic microscopic calculations including all previously 
mentioned effects are required. However, this is unconceivable in finite nuclei at this stage. 
As a first step, performing systematic HFB calculations with our effective interaction offers 
a unique opportunity to understand in detail the contribution of the bare force to pairing 
in nuclei. 

Let us now discuss the range of the effective pairing force. Indeed, it would be nice to 
have an interaction allowing for large scale microscopic calculations of nuclear masses [35, 
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37]. Performing sclf-consistcnt mean- field or beyond-mean- field calculations, the size of 
the single-particle basis necessary to get converged results should be as small as possible 
to make such large scale calculations tractable. The required size is directly related to 
the range of the force in the pairing channel. For instance the shorter range of the Gogny 
force (0.7 /m) makes the convergence of the quantities related to pairing quite slow. One 
needs to include states up to 100 MeV in the quasi-particle spectrum to do so [111]. It 
can be seen from Fig. 3 that our effective interaction is softer than the Gogny force. The 
range to be compared with the 0.7 /m of Gogny is in the present case V^a 1/m. This 
translates into a kinetic energy of about 70 MeV which, once the depth of the single- 
particle potential has been subtracted 50 MeV), gives a value of the order of 20 MeV. 
Taking the effect of the cut-off {pj + pi) in the gap equation into account, it should be 
sufficient to include single-particle states up to 10 — 20 MeV of positive energy in the 
canonical basis to get converged HFB calculations. Thus, this pairing interaction should 
be simple and soft enough to perform large scale calculations of nuclear masses, avoiding 
at the same time the problem related to the phenomenological choice of the cut-off dealing 
with the zero range of the DDDI [37] . Checking this statement is the aim of a forthcoming 
publication [117]. Note that the induced interaction is expected to be long ranged [52, 55]. 

Once the softness of the pairing interaction has been established, it will be important 
to understand the significance of the range and the non-locality of the force and whether 
it needs to be treated explicitly in nuclei. It amounts to studying the influence of the 
non-locality of the pairing field. Being able to use the microscopic non-local finite-range 
interaction and a gap-equivalent local zero-range force within a single scheme offers a 
unique opportunity to answer such a question [117]. 

5 Conclusions 

In this paper, we have proposed a microscopic effective interaction to describe pairing in 
the ^Sq channel. The new features of the interaction are numerous. 

It possesses a clear link to the bare force. The effective interaction provides the same 
pairing properties as the bare force at the mean-field level, as required from many-body 
theories. The gap at the Fermi energy obtained in infinite matter from the realistic 
AV18 interaction is perfectly reproduced by the new force. Going into further detail, 
the momentum dependence of the gap is also very well described in the density range of 
interest. 

The effective force is finite ranged, non local, total-momentum dependent and density 
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dependent. While the effective interaction is almost constant for densities ranging from 
saturation to typical surface densities, it is strongly enhanced at very low density. The 
isoscalar and isovector density dependences of the pairing force are also obtained through 
the ab-initio derivation. While phenomenological density-dependent pairing interactions 
used so far depend on the total density, the one derived here depends on the density of 
the interacting nucleons only (i.e. protons or neutrons since the present work deals with 
like-particle pairing only). 

This effective pairing force is defined by recasting the gap equation written in terms of 
the bare force into a fully equivalent pairing problem. Through this rewriting procedure, 
the matrix elements of the effective force are provided with a natural cut-off 2/9^, where 
Pm is a BCS occupation number. This makes the definition of zero-range approximations 
meaningful and no ad-hoc cut-off has to be additionally chosen and optimized. Performing 
such a zero-range approximation and asking for identical pairing gaps at the Fermi surface 
in infinite matter, the appropriate density-dependence of the zero-range force is obtained. 
This procedure, free of any phenomenological cut-off allows us to disentangle the roles 
played by the range and the density dependence of the pairing interaction. Surprisingly, 
the enhancement of the force at typical surface densities as compared to the saturation 
density is much more pronounced than for the finite range vertex. This result shows 
unambiguously that the surface character of usual zero-range forces is, to a large extent, 
a way of taking care of the range of the interaction. Precisely, the zero-range vertex 
is predicted to have a behavior between surface and volume. It also undergoes a large 
increase of its intensity at very low density. 

The last essential feature of the force can be identified when going from infinite matter 
to finite nuclei. Indeed, dealing with a finite-range and non-local interaction is far from 
being trivial when solving the HFB equations in coordinate space. However, the particular 
analytical structure of our effective force makes such calculations possible within the so- 
called two-basis method. In fact, the corresponding computational cost is of the same 
order as for zero-range forces. Performing exploratory calculations in finite nuclei with 
this new interaction is the aim of a forthcoming publication [117]. 
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